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ABSTRACT 



Humidity spectra have been measured with the Lyman - 

alpha humidiometer , together with mean profiles of virtual 

potential temperature, wind, and humidity in the open 

ocean environment. Empirically derived expressions des- 

2 

cribing the temperature structure parameter, , were 

extended by similarity arguments to the humidity-structure 
2 2 

parameter, C , and C was related to the stability para- 
q q 

meter Ri. Using the above measured parameters, vertical 

humidity flux was computed in two different manners , and 

a comparison was made. In general, there was little 

2 

correlation between the spectrally analyzed values 

2 

and Ri. Results for can essentially be regarded as 

a function of z and humidity gradient. Non-dimensional 
2 

0^ results were generally an order of magnitude smaller 
than expected. No correlation between the two methods of 
calculation of humidity flux was found. 
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I. INTRODUCTION 



In recent years the turbulent structure of the atmos- 
pheric boundary layer has become very important to both 
electro -magnetic and optical wave propagation studies. 

Small scale temperature and humidity fluctuations are 
important in determining such phenomena as the intensity, 
fluctuation, beam wander, beam spread, and absorption of 
a transmitted wave. It is also recognized that these 
fluctuations are important in the understanding of electro- 
magnetic waves. 

Initial experimental efforts to verify turbulence theory 
predictions were conducted over land. Measurements of wind 
speed, temperature and humidity fluctuations in those in- 
vestigations were taken from a stable platform with instru- 
mentation not encountering as severe conditions as were these. 
The Navy, however, is interested in the prediction of these 
parameters in the marine environment. Experiments at the 
Naval Postgraduate School have been providing the basis for 
predictions of the properties of the atmosphere for rele- 
vant wavelengths. The turbulent fluctuation of temperature 
is that parameter which is important in optical wave propa- 
gation. Numerous studies of the turbulent temperature 
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fluctuations have yielded significant results. In this 
analysis, another facet of the atmospheric boundary layer 
turbulent regime is considered, namely humidity. 

It was the purpose of this thesis to examine the turbu- 
lence properties of water vapor in the atmosphere. It was 
examined in a fashion similar to that used for temperature, 
since both parameters are passive scalar contaminants in 

the atmosphere. An analysis of the humidity-structure 

2 2 
function, C , was made, and an attempt to relate C to a 

q q 

stability parameter, the Richardson number was made. Like- 
wise, the turbulent transfer of the vertical humidity flux 
(w'q') was examined. 

This study was based on data obtained during April 1976 
under open ocean conditions aboard the R/V Acania. The 
data were examined on the basis of expressions used 
previously by Wyngaard et al. (1971) and Friehe et al, 
(1975) . 
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II. THEORETICAL BACKGROUND 



A 0 GENERAL 

Based on analyses of Corrsin (1951) the variance spectra 
of a passive scalar property, such as humidity, at large 
Reynolds numbers should have an inertial subrange of the 
form 

s 00 - $ X c‘ 1/3 k' 5/3 , (1) 

q q 

where is an empirical constant evaluated to be 0.25. 

The temperature spectrum in terms of the structure function 

2 

parameter C^ , has the form 

S T (k) = 0.25 C T 2 k‘ 5/3 . (2) 

Similarity arguments lead to a parallel expression for the 

2 

humidity structure function parameter, C^ 

2 -5/3 

S (k) = 0.25 C k ' . (3) 

q q 

Equation (3) defines the one -dimensional humidity spectrum 
which by definition is the Fourier transform of the corre- 
lation function with separation r in the streamwise 
direction. 
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B. CALCULATION OF THE HUMIDITY -STRUCTURE PARAMETER 

2 

Equation (3) can be solved for C as 

q 

2 5/3 

C = 4 S (k) (4) 

q q - 

where k is the streaniwise component of wave number. Since 

the humidity fluctuations are measured at a fixed point in 

the flow, the resultant spectra are realized at a temporal 

2 

frequency, f. To obtain C^ , the temporal (f) and space 
(k) scales are assumed to be related by Taylor's frozen- 
field hypothesis, 

k = 2 7T f /U (5) 



where U is the mean wind at the measurement level. "Frozen" 
implies that the turbulence pattern remains unchanged as it 
sweeps past the probe. This yields 



2 

C 

q 



- 4 S (k) (is-) 5/3 
q u 




( 6 ) 



This is the expression that was used to compute the humidity- 
structure parameter from humidity variance spectra, S^(k) . 



C. STABILITY CONSIDERATIONS 

Present boundary layer turbulence theory has its basis 
on empirical results obtained by Monin and Obukhov (1954). 
The Monin-Obukhov theory defines a scaling length, L, 
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proportional to the level where mechanical and thermal 
production of turbulent kinetic energy are equal. The 
Monin-Obukhov length is defined as 



L 



TyU * 2 

q * T. 

1 ^ V 



( 7 ) 



where T* v = T* + 0.61 T, q^ = -w'q'/U*, and K - 0.35 
(von Kantian' s constant). The ratio of the height of the 
measurement to the Monin-Obukhov length, z/L, serves as a 
stability index. 

Recent observational experiments by Businger et al 
(1971), as well as others, have provided a relationship 
between the Richardson number, Ri, 



Ri 



- g( ^Qv/ dz) 

du/ d z) z 



( 8 ) 



and the Monin-Obukhov length, L, where 9v is the virtual 
potential temperature. Figure (1) from Businger et al 
(1971) illustrates this determination quite well. The 
latter results substantiated expressions for the relation- 
ships between z/L and Ri proposed by Dyer and Dicks (1970) 
for unstable conditions. These expressions are, respectively, 
z/L = Ri (9) 



z/L 



Ri 

1- a Ri 



( 10 ) 
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Here d is an empirically derived constant equal to 0.5. 

It is important to note that in these expressions z/L 
approaches « as the Richardson Number approaches a critical 
value of 0.21. This suggests that as stability increases, 
the flow becomes essentially non-turbulent and the effect 
of mechanical turbulence becomes negligible. 




F igure 1 . 

The dependence of Richardson Number on stability. 
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Applying the definition of Eqn (4) to the general 



expression of Eqn (1) , it is seen that 



C 




(ID 



q 



q 



2 



0^ as defined by Eqn (11) enables indirect estimates of 
2 

C to be obtained based on mean conditions. This is because 

q 

e and X are easily related to boundary fluxes and profiles, 

if steady and horizontally homogeneous conditions exist. 

2 

Expressions which relate to mean properties of the bound- 
ary layer, such as Ri or z/L, are desirable because the small 
scale measurements are impractical to obtain in most opera- 
tional or tactical situations. 

2 

Expressions relating to € and X , and similarity 
theory predictions for the dependence of € and X on 



momentum and humidity fluxes were obtained from similarity 
arguments by Wyngaard et al (1971): 



4 - 



€ 



f ]_( z/L) 



( 12 ) 





(13) 




and f 2 (z/L) are empirically determined functions. 
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Direct substitution of Equations (12) and (13) into 



Equation (1) yields an expression of the form 
C q 2 = q* 2 z" 2/3 f 3 (z/L) . 



(14) 



f^(z/L) is defined by the combination of f^(z/L) and , 

and is given by 



f 3 (z/L) 



4.9 [1 - 7(z/L)] " 2/3 ; z/L < 0 

(15) 

4.9 [1+2. 8 (z/L)] ; z/L > 0 



Furthermore, since z/L and Ri are functionally related, 
as described previously, a parallel dependence of Ri can 
be obtained, which is (Wyngaard et al) 

C q 2 = z^ 3 ( dq/ Sz) 2 f^(Ri) (16) 

This final expression provides a desired dependence of 

2 _ 

C q on more readily measured parameters (z, dq/ d z) , and 

Ri) . 

Results for the non-dimensional temperature structure 
function parameter derived over open ocean conditions by 
Hughes (1976) and its comparison with the Wyngaard et al 
prediction curve are shown in Figure 2. 
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(Z P/9P) C/ z/ ^0 
z 'b ' z 
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Results of temperature-structure function parameter vs. Richardson Number (Hughes, 1976) 



D. CALCULATION OF VERTICAL HUMIDITY FLUX 



It is possible to determine the vertical humidity flux 



(w'q') from measurements of wind, C^ , and vertical humidity 
gradient. 

1. The Spectral Method 

The spectral method yields the following relation- 
ship for the viscous diffusion of humidity, X , 



X = (w'q’) d q/ d z 

q 



(17) 



Here, the humidity variance balance is determined for steady 
state, horizontally homogeneous conditions, and q is the 
specific humidity. 

Combining Equations (11) and (17) and solving for 
(w'q') yields 



..l-l 

w q 



2 1/3 

C q 6 

0.25 dq/dz 



(18) 



Eqn (18) needs to be simplified to use with the measured 
data. This simplification will be described in the next 
section. 

2. Direct Profile Method 

The acknowledged expression for the gradient of 
humidity pertinent to a steady-state, horizontally homo- 
geneous atmosphere can be written as: 
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d q / 3 z - q* / k z 



(19) 



Applying the definition of the scaling humidity, namely 



q* 




( 20 ) 



to Eqn (19) results in 



d q/ d z 



-w'q ' 

K z IL. 



( 21 ) 



When expressed in finite difference form, a simple approxi- 
mation to dq/^z is 



A3. 

d Z 



~ Aa = q.(. z 2) -q( . z i ) 

z 3 d z z 31n (z 2 /z 1 ) 



( 22 ) 



where the computed number applies to the height z 3 . Apply- 
ing this finite difference approximation to Eqn (21) and 
solving for (-w'q') yields 

= *u*(g(z 2 >-q( z i)) (23) 

In (z 2 /z 1 ) 

This is the expression that was used to compute vertical 
humidity flux from the vertical profile of mean specific 
humidity. 
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III. THE EXPERIMENT 



A. THE PLATFORM AND LOCATION 

Observations were made aboard the R/V ACANIA anchored 
off Monterey, California in Monterey Bay. Measurements were 
made at multiple levels on two masts spatially separated on 
the forward deck of the ship. An illustration of the sensor 
locations of these masts is shown in Figure (3) . The vane 
and probe arrangement appears in Figure (4) . The vane main- 
tained a one -dimensional profile of the wind during 
measurements . 

The location of Monterey Bay provided an ideal site for 
the experimental program as shown in Figure (5). Open ocean 
differs from land in the effects of wave action on turbu- 
lence, in the nature of the aerosols and fog, and in the 
humidity fluctuations. These conditions can, of course, be 
best obtained far at sea. The cost of such activities makes 
it desirable to work near land. Pt. Pinos and Monterey Bay 
provide very nearly the ideal situation. 

Pt. Pinos projects northwest from the mainland toward 
the prevailing northwest wind. Even under storm wind condi- 
tions, the wind comes from the southwest, still bringing 
sea air toward shore. The ocean depth surrounding Pt. Pinos 
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Figure 3. 

Mounting arrangements aboard the ACANIA 
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Vane and Probe Arrangement 






F igure 5 . 

Experimental site in Monterey Bay 
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is shallow enough that the ACANIA can anchor at any range 
up to about 15 km to the north and northeast. 

B. INSTRUMENTATION 

1. Wind-Temperature Measurement Systems 

The mean wind measurements were made with a 
Thomthwaite Associates cup anemometer wind profile register 
system, model number 104; In operation the shaft of a three 
cup anemometer unit serves as the shutter between a light 
source and photocell for each revolution. The cups are 
plastic cones reinforced with aluminum frames. They are 
attached to the rotating shaft by stainless steel tubing 
spaced at 120 degree intervals about the shaft as shown in 
Figure (6) . The three cup assembly sets along with the 
other sensors were positioned at four levels on the bow 
with electrical leads to the after deck house laboratory. 

The sets have the characteristics of low starting speeds 
with a small amount of internal friction which aids in 
checking inertial overshoot. 

Temperature sensitive quartz crystal probes, (Hewlett 
Packard model HP-2850) were used to measure mean temperatures 
at each level. RF signals from the crystal probes and from 
a reference oscillator were mixed in the HP-2801A readout 
unit to produce a beat frequency whose signature can be 
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F igure 6 . 

C. W. Thomthwaite Anemometer Cups 
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analyzed to within 0.001 degrees centigrade per hertz. Each 
sensor simultaneously received pre-experiment calibration 
against a platinum resistance wire thermometer in a tempera- 
ture controlled circulating water bath over the expected 
temperature range. The accuracy in achieving a 0.005 degree 
centigrade correction factor was a constant for each probe. 

A 3.7 meter flexible coaxial cable is permanently attached 
to the sensor head and the mast mounted probes are housed 
in an aspirated shelter as depicted in Figure (7). Tempera- 
ture values were automatically recorded on a printer tape. 

The temperature fluctuations were measured using a 
bridge developed by personnel at GTE Sylvania, the GTE 
Sylvania Model 140. The system was slightly modified for 
use in this study, and is described by Karch (1976). 

2 . Humidity Measurement Systems 

To determine their accuracy and reproducibility, 
detailed examination in laboratory conditions has been made 
of the calibration characteristics of two humidity sensors, 
the Hygrodynamics Digital Hygrosensor and the Electro- 
Magnetic Research Corporation Lyman-alpha Sensor. These 
calibrations were based on ambient humidity variations 
(from 24%, RH to 507, RH over a one month period) , as measured 
with a Bendix Psychron Wet-Dry Bulb Psychrometer , and con- 
trolled humidity (from 40% RH to 99%, RH) as measured by the 
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C.C. Breidert Company 
Air-X-Houster Type 6L 



Dunmore-type Lithium 
Chloride Sensor 




Quartz Thermometer Probe 
Figure 7. 

Photographs of aspirated shelter, hygrometer 
and thermometer probe 
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wet-dry bulb thermocouple in a 2 x 3 x 3 ft. fog chamber 
adapted for these experiments. The air circulated through 
the chamber can be pre-humidif ied in a column of wet glass 
beads giving nominal humidity of about 80%. Higher humidi- 
ties (up to 99% RH) can be obtained by introducing varying 
amounts of water in the form of a fine mist into the chamber. 
In general, the humidity in the chamber can be stabilized for 
periods of several hours and is uniform throughout the 
chamber to within 1%. RH. 

a. Hygrodynamics Digital Hygrosensor 

The hygrosensor is a Dunmore-type lithium 
chloride sensor whose resistance varies in proportion to 
the relative humidity to which it is exposed. It is a s low- 
response sensor and is generally utilized in mean system 
measurements. The standard procedure of calibrating these 
sensors in a small chamber with a saturated salt solution 
of known vapor pressure was not satisfactory for a program 
of calibration of eight sensors. This dissatisfaction was 
primarily due to a lack of confidence in the ability to 
reproduce humidity conditions in consecutive sensor calibra- 
tions. Attempts to calibrate all eight sensors simultaneously 
in a large volume chamber were complicated by temperature 
drifts that made it difficult to assure equilibrium satura- 
tions necessary to the accuracy of the technique. 
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As a result of the above problems it was decided 
to use the humidity controlled chamber previously described. 
The quantitative results of their calibration are given in 
Table I; in general, the eight sensor average agrees with 
the psychrometer standard to -0.4 + 2.9% RH for measurements 
taken with equilibrium times of order one hour. The indi- 
vidual total humidity accuracy of the sensor is of less 
importance that the relative consistency within the group 
of sensors, since we use them in a four-level system to 
determine the humidity gradients. In this respect, they can 
be calibrated to + 1% RH, subject to assurances of sufficient 
time for response. 

The suitability of the hygrosensors for ship- 
board measurement of humidity gradients is reduced by two 
inherent properties of the device. An aspirated, teflon- 
coated sensor has a time constant of about 30 minutes, but 
this varies from sensor to sensor. Therefore, the humidity 
difference read between two sensors will be unreliable if 
the humidity is changing faster than about 10 % RH per half- 
hour. A more serious deficiency appears if the sensors are 
exposed to humidities about 95% RH. At those high humidi- 
ties, the sensors become sluggish and exposure to 99% RH 
may cause the sensors to be useless for as long as several 
days . 
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TABLE I 



Hygrodynamics Digital II Calibration 



Wet- 

Dry 


Ave 


Deviation of 


Bulb 

%RH 


Hygro 

%RH 


1 


2 


3 


33. 


3.16 


1.4 


-1.4 


-0.6 


36. 


37.2 


-0.5 


-0.8 


-0.3 


44. 


41.2 


-0.8 


-0.4 


-0.3 


NA 


41.5 


-0.5 


-0.4 


-0.3 


NA 


42.6 


-1.2 


-0.5 


-0.5 


NA 


65.8 


-0.2 


-1.4 


-2.5 


74.8 


73.6 


0.2 


0.2 


-1.3 


88.3 


85.0 


-0.7 


0.9 


-0.4 


82.6 


86.7 


-1.4 


-0.6 


-0.5 


82.6 


87.0 


-1.0 


0.5 


-0.7 


85.0 


88.0 


-0.4 


1.1 


-0.1 


93.9 


89.6 


0.1 


0.8 


-0.2 


94.4 


92.3 


-0.6 


1.1 


-0.5 


94.4 


93.6 


-0.8 


1.0 


-0.8 


Ave 


Deviat 


-0.5 


-0.0 


-0.6 






0.7 


0.9 


0.6 



Sensor from Ave Hygro % RH 



4 


5 


6 


7 


• 8 


1.5 


-0.2 


-1.6 


-0.9 


1.8 


1.6 


0.5 


-1.4 


-0.1 


1.2 


0.6 


0.9 


0.4 


-0.1 


0.3 


0.5 


0.9 


-0.4 


-0.1 


-0.1 


-0.9 


1.8 


-0.4 


-0.2 


0.1 


0.8 


1.8 


-0.3 


1.3 


0.3 


0.5 


-0.2 


0.4 


-0.4 


0.8 


0.1 


-0.6 


0.0 


0.4 


0.2 


0.1 


-0.2 


0.3 


0.7 


0.3 


0.0 


0.3 


0.7 


0.0 


0.2 


0.1 


-1.8 


1.1 


-0.6 


-0.1 


0.3 


-2.2 


2.1 


-1.5 


0.4 


-0.3 


-0.4 


1.2 


-0.5 


0.4 


-0.3 


0.4 


-0.2 


0.8 


-0.4 


0.3 


0.1 


0.1 


-0.1 


0.4 


0.7 


1.1 


1.0 


0.7 


0.6 



The average deviation of the wet -dry bulb from the Ave Hygro 
reading is 0.4% RH with a standard deviation of 2.9% RH. 
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b. Electro-Magnetic Research Corporation 
Lyman-alpha Sensor 

The Lyman-alpha sensor measures humidity as a 

function of the absorption of ultraviolet light by the Model 
o 

1215A Lyman-alpha transition of hydrogen in water vapor. It 
consists simply of a UV source tube and detector, separated 
by an absorbing air gap. 

The detector voltage is proportional to the 
total light transmitted and is related to the absolute 
humidity, q (in mb) , by the equation 



V = V Q exp ( - A q ) 



(24) 



The Lyman-alpha is a fast response device and is specifically 
used for measurement of humidity fluctuations, q'. Eqn (24) 
can be manipulated to relate q' to voltage fluctuations V', 
giving 



1_ VI 
A V 



(25) 



The quantity V in Eqn (24) is the sensor output voltage in 
o 

dry air which can be found by placing the detector in de- 
humidified air. 

Because water soluble windows (LiF or MgF^) are 
required for transmission in the ultraviolet, V q is a highly 
variable quantity, a property which makes the Lyman-alpha a 
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poor device for measurement of total humidity, q. However, 
Eqn (25) shows that the fluctuating part, q', depends only 
on A , which is a property of water vapor and the gap setting 
only and is independent of the window condition so that q' 
can be reliably measured. The Lyman-alpha has been cali- 
brated in dry air, ambient air, and in a humidity chamber 

and it was found that A = .225 mb“^ for a source current = 

50 /ia and a 1 cm air gap. The calibration curve appears 

in Figure (8) . 

The feasibility of the Lyman-alpha sensor for 
shipboard humidity fluctuation measurements is subject to 
considerations of the survivability of the windows in the 
ocean environment of high humidity, sea spray, and rain. 

In the laboratory humidity chamber, it was observed that 
V Q would decrease by a factor of four after a three-day 
exposure to 857, humidity. Consequently, it will be neces- 
sary to provide better physical protection for the sensors 
during routine operations. 

3 . Data Recording 

All mean data recorded for this study was logged 
using an NPS developed micro-processor based data acquisi- 
tion system, the MIDAS (Mic reprogrammable Integrated Data 
Acquisition System) . It is fully described by Atkinson 
(1976). In this study a ten minute averaging period was 
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V (Volts) 




Figure 8. 

Calibration graph of Lyman-alpha humidiometer 
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used to define the mean data. Temperature, relative humidity, 
and mean wind data was punched on paper tape by the MIDAS 
unit; the data was later transferred to computer cards for 
further analysis. 

Voltage signals corresponding to measured fluctuat- 
ing humidity, q', were recorded on magnetic tape using a 
Sanborn Model 3950 fourteen channel tape recorder. In addi- 
tion, the average voltage output of the Lyman-alpha sensor 
was recorded on a brush type chart recorder. The latter 
was necessary to define V in the voltage to humidity con- 
version expression, Eqn (25) . 

C. ANALYSIS PROCEDURES 

Data analysis was carried out for data collected during 
a four day cruise in April 1976. Mean data was obtained for 
each of the four levels during the entire cruise. However, 
only two Lyman-alpha humidiometers were available so q' data 
was collected at only two levels. These two instruments 
were placed at levels 6.6 and 13.9 meters above mean water 
level. Due to the failure of a source tube in the instru- 
ment at the upper level, q' data was available from only 
the 6.6 meter level. Table II summarizes periods examined 
and data available for each period. 



37 



TABLE II 



Summary of Data Periods 



Date 


Time 


Number of 
Levels 

_ q/q ' 


Number of 
Profiles /Spectra 

_ q/q ' 


1976 

27 April 


1446-1930 


4/1 


12/6 


28 April 


1320-1939 


4/1 


33/24 


29 April 


1410-2240 


4/1 


42/24 


30 April 


1033-1138 


4/1 


21/13 


All data 


were first 


edited for gross 


errors or incon- 



sistencies due to equipment malfunctions. The criterion at 

this point for retaining or rejecting data periods depended 
2 

on whether Cq values could be computed from humidity spec- 
tra for the time periods involved. 

Lyman-alpha data were analyzed on the basis of the strip 
charts to give an average voltage output at ten minute inter- 
vals. Relative humidity data from the MIDAS output was 
converted to specific humidity q, by use of the governing 
meteorlogical equations for moist air. Virtual potential 
temperatures, 9v, were calculated using IBM 360 routines. 

The values of q and 0v were then plotted versus the logarithm 
of height. Best fit lines were drawn to the points from 
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which slopes were picked off and applied to the expressions 
developed by Wyngaard et al. 

1 . Profile Analyses 

Virtual potential temperature, 0v, and specific 
humidity, q”, values were plotted on semi-logarithmic paper, 
since 0v, and q" are parameters which vary logarithmically 
with height. A best fit straight line was then drawn to 
the data points. The procedures, of course, were subjective, 
and in many instances several different slopes could be 
obtained from just one graph. Hence a criterion used was 
not to give single data points too much weight in determining 
the best fit line. Consequently, the line drawn represented 
a most probable position between data points as illustrated 
in Figure (9) . 

Difficulty arose whenever anomalous points appeared 
in a graph as in Figure (10) where the specific humidity at 
the second level was obviously inconsistent with the other 
three levels. In a case such as this, the inconsistent 

point was merely ignored. 

2 

2. C q Analysis Procedures 

Analog spectral analyses procedures were applied to 
the humidity fluctuation data to obtain variance spectrum. 
Appropriate 21 minute segments of data recorded on magnetic 
tape were recorded into an EMR-Schlumberger Model 1510 



39 



Z(m) Z(m) 





©v(°Cl 



Figure 9. 

Typical profiles for 

a) specific humidity - 30 April 76 (1210) 

b) virtual potential temperature - 29 April 76 
(1052) 
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9 (9/kg) 



Figure 10. 

Profile of 27 April 76 (1817) showing an 
anomalous specific humidity reading at the 
6.6 meter level. 
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Digital Spectrum Analyzer. The procedures required to con- 
vert spectral values, obtained with the digital spectrum 

2 

analyzer, to engineering units and then to obtain values 
from appropriate spectra are described in the following 
subsection. A parallel discussion for temperature spectra 
appears in Hughes (1976). 

a. Scaling Spectral Plots 

A necessary procedure was to scale the spectral 
plots to relate RMS input voltages to power spectral densi- 
ties (PSD) ; variance per unit frequency. To obtain power 
spectral density levels, corresponding to RMS voltage inputs, 
calibrated scale charts had to be constructed. 

For purposes of the x-y plot format of the 
analyzer output, the RMS voltages were converted to y=log-^Q 
(voltage) units and a graduated scale was constructed so 
that the logarithm of volts RMS could be interpolated from 
spectral plots. The value of the vertical scale (y) was 
adjusted for each spectrum as a function of both input gain 

and spectral gain. These values were then converted to PSD 

2 

levels for use in calculating values. 

The relation between Volts RMS to PSD units 



used was 

S(f) {PSD units} 



(cal, level VRMS) 2 „ no y ) 2 
1 . 5 Bandwidth 



( 26 ) 
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where Bandwidth = ^ r — ^ 

No. Channels 



/ _ 256 Hz _ 
V 256 



1 Hz) 



and cal. level Vrms = voltage at y=0 

( = 1 Vrms for 3.16 V input 
setting) . 



Amplitude scaling calibrations were accomplished 
using externally generated "white noise" signals of 1 Volt 
RMS with a frequency range from 0.1 Hz to 1000 Hz (giving a 
PSD of 10~ 3 V 2 /Hz). Setting a 3.16 V (lOdB) input on the EMR 
1510 Digital Spectrum analyzer insures that 1 Volt corre- 
sponds to y=0. An example of such a calibration plot for 
a frequency range of 1-200 Hz (the frequency range of inter- 
est) is shown in Figure (11) . 

b. Computation of Turbulence Parameters from 
Scaled Spectra 

2 

Values of the turbulence parameter C^ were 
obtained from the humidity variance spectra on the basis 
of the formula for the inertial subrange in wave number 
space S(k) , Eqn (1), which predicts a -5/3 slope for the 
spectra when plotted in a log-log format. Figures (12) , 
(13), and (14) are typical spectra considered in the analy- 
ses. Humidity spectra often exhibited slopes slightly 
different than the expected -5/3, except in the frequency 
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Figure 11. 

Calibration plot for spectrum analyzer 
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Figure 12. 

Typical humidity spectra showing -5/3 
slope in the range from 1 Hz to 10 Hz 
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Hz 



Figure 13. 

Typical humidity spectra showing -5/3 
slope in the range from 1 Hz to 10 Hz 
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.5 




Figure 14. 

Typical humidity spectra showing -5/3 
slope in the range from 1 Hz to 10 Hz 
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range of 1 Hz to 10 Hz. This feature of humidity spectra 
has been observed by others, namely Friehe (1975) . A uni- 
versally accepted explanation has not been postulated to 
account for this phenomena. Figure (15) is an example of 
a spectra which does not exhibit a -5/3 slope in any fre- 
quency range. Approximately 107, of spectra studies fell 
into this category. 

Assuming -5/3 slopes for the variance spectra, 

the intercept of the best fit -5/3 slope with the 1 Hz 

frequency line was the spectral density denoted (PSD) value 

2 

used to compute . The measured PSD value was converted 
to engineering spectral density units by the relation 



S (f) = C T 2 * PSD 
q v y L- a 

» (g/kg/Volt) 2 • Volt^/Hz = (g/kg) 2 /Hz 

(27) 

where was the calibration factor for the Lyman-alpha 

humidiometer . 

Since the humidity fluctuations were measured 

at a fixed point in the flow, the resultant spectral values 

are defined at "temporal" frequencies denoted as S^(f) in 

2 

Eqn (28). In order to obtain from Eqn (6) temporal (f) 

and space (k) scales had to be related. This was accomplished 
by using Taylor's "frozen turbulence" hypothesis discussed 



48 




Figure 15. 

Humidity spectra lacking a distinct -5/3 slope 
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in section II. B. The following equation relates temporal 
to wavenumber variance spectral values, 

fs q (f) - ks q (k) = c 2 C q V 2/3 (28) 

where Sq(f) is the spectral density value with units of 
(g/kg) ^ /Hz . 

Eqn (28) leads to the following relationship 
between and S^(f) , 

2/3 

c q 2 = !s— fS q (f) (29) 

where, C 2 = 0.25, an empirical constant, and k = 2 7rf/U. 

_ 2 

Hence, with measured values of f, S^(f), and U, Cq was then 

determined from Eqn (29) for the level of interest, z. 

3. Procedures to Calculate Vertical Humidity Flux 

Vertical humidity flux (w'q') was calculated using 

two different methods ; one utilized the spectrally derived 

2 

Cq values, and the other utilized the profile of mean 
humidity. The results were compared on a linear plot to 
determine if any relationship could be observed, 
a. The Spectral Method 

This method of calculating (w'q 1 ) was based on 

2 

Cq values derived from variance spectra of the Lyman-alpha 
humidiometer using Eqn (18). In this calculation the values 
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of € which were used were derived from wind fluctuation 
data by Karch (1976) . 

b. The Profile Method 

This method utilized the vertical gradient of 
mean humidity and the friction velocity, ILl, to calculate 
(w'q') from Eqn (23). The profiles of q" were evaluated 
and values of q were read off at two different levels from 
the best fit line drawn to all but erroneous points. 

The values of U* used in these computations 
were determined from the € values obtained by Karch (1976). 
The final relation used was 

U* - ( € k z) 1/3 , (30) 

which appears to give the most consistent estimates for Uj. 
in the marine boundary layer. 
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IV. RESULTS 



Observed humidity index- structure parameter results, 
made non-dimensional according to the Wyngaard et al predic- 
tion, Eqn (16), appears in Figure (16). In the figure 
individual data points appear as dots , and averages over 
Ri intervals of 0.25 appear as dots within a larger circle. 
The error bars are standard deviations from the mean within 
each interval, while the number at the top of the error 
bars is the number of observations defining the mean value. 
For both the stable (+Ri) and unstable (-Ri) stratifica- 
tion cases, there appears to be little agreement and no 
definite trends. 

Scatter in the observed results in Figure (16) can be 

o 

attributed to scatter in both the measured values as 

well as dq/ dz values. Deviation of humidity spectra from 

2 

a -5/3 slope also caused uncertainty in estimates. 

The values for d qV d z are important because the humidity 

gradients were, in most cases, small, and enter into the 

2 

normalization of as a squared value. 

2 

Significantly, non-dimensional results obtained in 

this study were generally an order of magnitude smaller 
2 

than Gj, results obtained by Wyngaard et al (1971) and 
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Figure 16 

Results of the spectrally derived humidity- structure 
parameter vs. the Richardson Number 
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and Hughes (1976). All sensors, calibrations, and govern- 
ing equations were rechecked. All calculations were redone. 
However, the final results remained smaller in value than 
expected. It was expected that results would fall in the 
same order of magnitude as C ^ since both temperature and 
humidity are scalars. One possible explanation is that the 
value of /? from Eqn (1) is not the same for both tempera- 
ture and humidity, as presently postulated. 

2 

If the scatter of the non-dimensional values is 

ignored, the results can be seen to be essentially constant. 

If one is willing to accept the error of an order of magni- 

_ 2 

tude, then given a gradient of q and a height z, can be 

computed according to the expression 

o 4/3 2 

C q Z = t 7 (dq/dz) Z (31) 

The comparison of the vertical humidity flux calculations 

n 

appear in Figure (17). (w'q') obtained from data are 

denoted (w'q')^_ a and vertical humidity flux obtained from 
profiles of dc[/ dz are denoted (w'q')p. There is no ob- 
served correlation between the two values. Perhaps new 
dimens ionalizing parameters need to be developed in order to 
be able to predict values of (w'q') accurately. 
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Figure 17. 

Results of the comparison between vertic al hu midity 
flux calculated by the profile method, (w'q')p, vs. 
the vertical humi dity flux calculated using 
C 2 values, (w'q') L _ a 
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Table III summarizes data used to construct Figure (16). 
Tables IV and V contain values used in the comparison shown 
in Figure (17) . Values in all three tables were computed 
at the 6.6 meter level. 
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TABLE III 



C 



2 



q 



and Ri Results 



Time 


U 


6q/5z 


C q 2 X lCf 3 


Ri 


27 April 
1446-1509 


7.70 


-0.043 


1.48 


-1.26 


1509-1520 


8.17 


-0.045 


0.88 


-1.07 


1520-1540 


7.84 


-0.048 


0.92 


-1.44 


1540-1600 


7.50 


-0.037 


1.08 


-1.27 


1620-1640 


6.25 


-0.045 


0.96 


-2.73 


1910-1930 


4.03 


-0.030 


1.72 


-0.90 


28 April 
1320-1337 


5.16 




2.60 




1337-1355 


5.30 


-0.026 


2.04 


-1.32 


1355-1413 


5.16 


-0.025 


1.72 


-0.93 


1415-1420 


5.06 


-0.028 


1.24 


-1.05 


1433-1451 


5.02 


-0.024 


1.96 


-1.32 


1452-1504 


4.39 


-0.024 


2.08 


-0.95 


1511-1517 


4.88 


-0.023 


2.00 


-1.26 


1526-1533 


4.50 


-0.031 


1.84 


-2.96 


1540-1551 


4.31 


-0.031 


2.60 


-3.20 


1556-1606 


4.40 


-0.024 


2.05 


-3.70 


1613-1621 


5.30 


-0.032 


1.67 


-1.42 


1633-1641 


4.97 


-0.039 


1.52 


-2.19 


1652-1701 


4.73 


-0.038 


1.27 


-0.58 


1720-1738 


4.31 


-0.044 


1.09 


0.19 


1739-1744 


4.03 


-0.051 


1.56 


-0.74 


1757-1807 


4.31 


-0.048 


1.00 


-0.58 


1807-1817 


4.31 


-0.053 


1.29 


-0.58 


1817-1827 


4.08 


-0.057 


1.36 


-0.52 


1827-1837 


4.03 


-0.052 


1.32 


0.24 


1837-1847 


3.91 


-0.052 


1.07 


0.29 


1847-1857 


3.61 


-0.053 • 


1.16 


0.29 


1909-1919 


2.80 


-0.045 


1.72 


1.68 


1919-1929 


2.50 


-0.033 


1.51 


7.01 


1929-1939 


2.50 


-0.028 


1.60 


7.89 


29 April 
1410-1420 


5.88 


-0.020 


2.36 


0.79 


1420-1430 


5.86 


-0.039 


3.36 


0.12 


1442-1452 


5.86 




6.40 




1500-1513 


6.43 




3.45 





Table III Continued 



Time 


U 


d q/ d z 


C q 2 x 10' 3 


Ri 


1521-1532 


6.00 




3.10 




1541-1553 


4.50 




3.77 




1600-1610 


5.86 




2.51 




1653-1703 


4.26 


-0.039 


2.36 


-6.15 


1703-1713 


4.92 


-0.024 


1.72 


-4.06 


1713-1723 


4.31 


-0.023 


1.68 


-10.94 


1723-1733 


4.22 


-0.014 


1.71 


-5.17 


1733-1743 


4.26 


-0.011 


2.60 


-4.82 


1830-1848 


2.60 




2.28 




1848-1906 


3.61 




2.41 




1916-1934 


2.71 




2.00 




1934-1952 


2.83 




2.16 




2000-2018 


2.95 








2018-2036 


3.36 








2036-2054 


3.76 




1.27 




2054-2114 


3.81 




1.25 




2125-2145 


3.66 




1.40 




2145-2157 


3.76 




1.60 




2203-2223 


3.76 




1.69 




2221-2240 


4.71 




1.60 




30 April 
1033-1043 


4.73 


-0.047 


2.04 


-1.56 


1043-1053 


4.45 


-0.043 


1.55 


-2.52 


1053-1103 


4.31 


-0.030 


2.44 


-3.91 


1108-1118 


4.40 


-0.043 


1.76 


-2.90 


1118-1128 


4.55 


-0.043 


1.75 


-2.52 


1128-1139 


4.78 


-0.033 


3.00 


0.37 


1148-1206 


5.53 


-0.031 


1.44 


0.74 


1206-1224 


5.53 


-0.036 


1.65 


0.62 


1224-1242 


6.05 


-0.054 


1.04 


0.0 


1242-1250 


4.40 


-0.041 


1.52 


0.56 


1250-1308 


3.26 


-0.034 


1.83 


0 . 0 . 


1310-1328 


3.41 


-0.040 


1.32 


-1.98 


1328-1346 


3.51 


-0.032 


0.83 


-0.96 
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TABLE IV 



w'a' Results. C ^ Method 



Time € (xlO 



28 Apr 76 
1337-1355 


2.90 


1355-1413 


3.67 


1757-1807 


4.55 


1807-1817 


4.60 


1817-1827 


4.66 


1827-1837 


4.76 


1837-1847 


4.95 



29 Apr 76 




1653-1703 


4.00 


1703-1713 


3.50 


1713-1723 


2.33 


1723-1733 


2.80 


1733-1743 


3.04 



30 Apr 76 




1033-1043 


0.97 


1224-1242 


2.09 


1242-1250 


2.09 


1250-1308 


1.37 



C q (xlO ) d q / d z 



2.04 


-0.026 


0.045 


1.72 


-0.025 


0.042 


1.00 


-0.048 


0.014 


1.28 


-0.053 


0.016 


1.36 


-0.057 


0.016 


1.32 


-0.052 


0.017 


1.16 


-0.053 


0.015 



2.36 


-0.039 


0.038 


1.72 


-0.024 


0.044 


1.68 


-0.023 


0.039 


1.72 


-0.014 


0.069 


2.60 


-0.011 


0.014 



2.04 


-0.047 


0.017 


1.04 


-0.054 


0.010 


1.52 


-0.041 


0.019 


1.84 


-0.034 


0.024 



59 



TABLE V 



w'q' Results, Profile Method 



Time 


q(z 2 ) -q(z 1 ) 


u* 


,, i _ i 

w q 


28 Apr 76 


1337-1355 


6.15-6.18 


0.19 


0.013 


1355-1413 


6.27-6.30 


0.20 


0.014 


1757-1807 


6.77-6.83 


0.22 


0.030 


1817-1827 


6.78-6.84 


0.22 


0.030 


1827-1837 


6.79-6.85 


0.22 


0.030 


1837-1847 


6.82-6.87 


0.23 


0.027 


29 Apr 76 


1653-1703 


6.30-6.35 


0.21 


0.024 


1703-1713 


6.41-6.43 


0.20 


0.009 


1713-1723 


6.39-6.42 


0.18 


0.013 


1723-1733 


6.33-6.35 


0.19 


0.009 


1733-1743 


6.40-6.42 


0.19 


0.009 


30 Apr 76 


1033-1043 


6.62-6.67 


0.13 


0.015 


1224-1242 


7.11-7.17 


0.17 


0.023 


1242-1250 


7.12-7.17 


0.17 


0.020 


1250-1308 


7.08-7.13 


0.15 


0.017 
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